NONLINEAR EVOLUTION PDES IN M+ x C'': EXISTENCE 
AND UNIQUENESS OF SOLUTIONS, ASYMPTOTIC AND 
BOREL SUMMABILITY PROPERTIES 
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Abstract. We consider a system of n-th order nonlinear quasilinear 
partial differential equations of the form 

Ut+P(^x)u + g(x,t,{c^,u}) =0; u(x,0) = u,(x) 

with u G C^, for t G (0, T) and large |x| in a poly-sector 5* in 

{d^ = dl\dl^2---'^'d ^'^'^ Ji + ■•• + id < n). The principal part of the 
constant coefficient n-th order differential operator V is subject to a 
cone condition. The nonlinearity g and the functions u/ and u satisfy 
analyticity and decay assumptions in 5*. 

The paper shows existence and uniqueness of the solution of this 
problem and finds its asymptotic behavior for large |x|. 

Under further regularity conditions on g and u/ which ensure the ex- 
istence of a formal asymptotic series solution for large |x| to the problem, 
we prove its Borel summability to the actual solution u. 

The structure of the nonlinearity and the complex plane setting pre- 
clude standard methods. We use a new approach, based on Borel- 
Laplace regularization and Ecalle acceleration techniques to control the 
equation. 

These results are instrumental in constructive analysis of singularity 
formation in nonlinear PDEs with prescribed initial data, an application 
referred to in the paper. 

In special cases motivated by applications we show how the method 
can be adapted to obtain short-time existence, uniqueness and asymp- 
totic behavior for small t, of sectorially analytic solutions, without size 
restriction on the space variable. 
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1. Introduction 

1.1. General considerations. There are relatively few general results on 
existence, uniqueness and regularity of solutions of partial differential equa- 
tions in the complex domain when the conditions of the classical Cauchy- 
Kowalewski (C-K) theorem are not met. The C-K theorem holds for first- 
order analytic systems (or those equivalent to them) with analytic non- 
characteristic data, and for these it guarantees local existence and unique- 
ness of analytic solutions. As is well known, its proof requires convergence 
of local power series expansions. Evolution equations with higher spatial 
derivatives do not satisfy the C-K assumptions and even when formal power 
series solutions exist their radius of convergence is zero. One of the goals 
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of this paper is to provide a theory for existence, uniqueness and regularity 
of solutions in such cases, in a relatively general setting. The theory also 
applies to classes of equations of higher order in time and sufficiently high 
order in space after reduction (by well known transformations, see e.g. [20!) 
to evolution systems. 

The present paper generalizes jHj to d dimensions and arbitrary order in 
the spatial variable, to r dimensional dependent variable, proves additional 
results about short term existence and shows Borel summability of formal 
solutions. A fortiori we obtain results on the asymptotic character of these 
solutions. (In Appendix H7.21 we briefly discuss the definition and properties 
of Borel summation.) 

Under assumptions to allow for formal expansions for large x, we show 
that series solutions are Borel summable to actual solutions of the PDE. For 
this purpose we make use of Ecalle acceleration techniques. In special cases 
we obtain existence and uniqueness results for t in a compact set and large 
enough x, and separately for small t and fewer restrictions on x. 

Properties of solutions of PDEs in the complex plane, apart from their 
intrinsic interest, are relevant for properties in the real domain, as initial 
singularities in C may give rise to blow-up at later times in the physical 
domain. Representation of solutions as Borel sums is instrumental in ex- 
tending techniques originally developed for ODEs [7j to find the location 
and type of singularities of solutions to nonlinear PDEs 

It is certainly difficult to give justice to the existing theory of nonlinear 
PDEs, and we mention a number of results in the literature relevant to the 
current paper. For certain classes of PDEs in the complex domain Sam- 
martino and Caffisch ^7j, proved the existence of nonlinear Prandtl 
boundary layer solutions for analytic initial data in a half-plane. This 
work involves inversion of the heat operator dt — dyy and uses the abstract 
Cauchy-Kowalewski theorem for the resulting integral equation. While their 
method is likely to be generalizable to certain higher-order partial differential 
equations, it appears unsuitable for problems where the highest derivative 
terms appear in a nonlinear manner. Such terms cannot be controlled by 
inversion of a linear operator and estimates of the kernel, as used in (|17j. 

CHI). 

The complex plane setting, as well as the type of nonlinearity allowed 
in our paper, do not allow for an adaptation of classical, Sobolev space 
based, techniques. This can be also seen in simple examples which show 
that existence fails outside the domain of validity of the results we obtain. 

Certainly, many evolution equations are amenable to our setting; to il- 
lustrate canonical form transformations and the general results we chose a 
third order equation with quartic nonlinearity arising in fluid dynamics. De- 
tailed singularity study (11; of solutions of this equation relies on the present 
analysis. 
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Our approach extends Borel transform regularization to a general class 
of nonlinear partial differential equations. A vast literature has emerged 
recently in Borel summability theory, starting with the fundamental con- 
tributions of Ecalle (see e.g. |13) ) whose consequences are far from being 
fully explored and it is impossible to give a quick account of the breadth 
of this field. See for example [Jj for more references. Yet, in the context 
of relatively general PDEs, very little is known. For small variables, Borel 
summability has been recently shown for the heat equation |16| , and gen- 
eralized to linear PDEs with constant coefficients by Balser [2j. One large 
space variable was considered by us in jHj, in special classes of higher order 
nonlinear PDEs. The methods in the present paper are different and apply, 
for large |x|, to a wide class of equations. 

1.2. Notation. We use the following conventions. For vectors in or 
multiindices we write 

d 

|u| = ^ \Ui\ 

and for multiindices we define 

k ;^ m if A;j > TTij for all i 

If a is a scalar we write x'^ = (x^, , x^). 

With p, x and j vectors of same dimension d, we define 

d 



and 



i=l 



We write 1 = (1, 1, .., 1) and more generally, if a is a scalar, we write a = al; 
thus x-*- = ni=i ^i- -^o^ d-dimensional vectors a and b we write 

fb r-bi r-b2 f-b^ 

• dp = / / ... / ■ dpidp2 ■ ■ ■ dpd 

The directional Laplace transform along the ray arg pi = ipi,i = l...d of F 
is given by 

(1) {C^F} (x) ^ / F(p)e-P-dp 

Jo 

where xe*^ will denote the vector with components XjC*^*. Convolution is 
defined as 



(2) 



if*9){p)-=[ /(s)5(p-s)ds 
Jo 



and Y\ denotes convolution product (see also ,6j). Whenever used as sum 
or product indices, I takes all integer values between 1 and m, i is between 
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1 and d, As a sum or product multiindex, |j| indicates all j with positive 
integer components subject to the constraint 1 < |j| < n. 

2. Problem statement and main results 

2.1. Setting and assumptions. Consider the initial value problem for a 
quasilinear system 

(3) ut + r{ci)u + g (x, t, {^xu}|j|<n) = 0; u(x, 0) = U7(x) 

In ©, 'P((9x)u collects the constant coefficient linear terms of the partial 
differential equation. 

Emphasizing quasilinearity, we rewrite the equation as 

(4) dtu + r{d^)u + g2,j (x, t, {^,u}|j|<„) diu 

J|=n 

= gi (x, t, {«9iu}|j|<„) ; u(x, 0) = U7(x) 

The restrictions on gi , g2 , and ui are simpler in a normalized form, more 
suitable for our analysis. By applying to (llj for all j with 1 < |j| < n — 1, 
we get an extended system of equations for f G C™, consisting in u and its 
spatial derivatives of order less than n, of the type (see Appendix for further 
details) : 
(5) 

atf + P(ax)f = j;'bq(x,t,f)J](^;,/,)'''+r(x,t) with f(x,0)=f7(x) 
where means the sum over the multiindices q with 

m 

(6) Yl Yl IjI^'j - ^ 

1=1 l<|j|<n 

The matrix V is assumed to be diagonalizable, and modulo simple changes 
of variables we assume it is presented in diagonal form, V = diagV j,j = 
l,...,m. In @, q = (qij), 1 < |j| < n, 1 < / < m is a vector of integers 
and Vj is an n-th order polynomial. We let Vnj be the principal part of 
Vj, i.e. the part that contains all monomials of (total) degree n. The 
inequality ® implies in particular that none of the g/j can exceed n and 
that the summation in ((S)) involves only finitely many terms. The fact that 
@ can always be ensured leads to important simplifications in the proofs. 
Let /?>/3o>0, (/><^,e>0 and 

(7) ^0,p;x = |x : I argXjl < ^ + (f>; \xi\ > p; i < d\ 



(8) 



'D^,P = 'D^,p;^ X [0,r] 
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Assumptions 1. (1) There is a (p ^ (O, ^) such that for a// p / with 
maxj I argpil < (j) we have 

(9) ^Vn;j{-P)>0 

(2) The functions bq(-,t,-) are analytic in ^^,po ^ {f '■ |f| < e}- 

write 

(10) bq(x,t;f) =^bq,k(x,t)f^ 

(3) For some constants ar > 1 independent ofT (see also ^7.1[ ), Ar{T) > 
0, aq > 0^ 



(11) 


sup 


|x"^r(x,t)| 


= Ar{T) < oo 


(12) 


sup 


|x*-f/(x,t)| 


= Af{T) < oo 


(13) 


sup 


Ix^'^bq^kl 


= Ab{T) < oo 



(4) T/ie analysis is interesting for n > 1, which we assume is the case. 
2.2. Existence and uniqueness for large |x|. 

Theorem 1. Under the AssumptionsU[ there is a unique solution f of |3|) 
satisfying the following properties in I',/,,po;x- (o-) f analytic and (b) |x-'-||f| 
bounded . Furthermore, this solution satisfies f = 0(x^"'') as x ^ oo in 
T^4>,p;^, for large p. 

Notes. 1. The proof of Theorem ^ is given in ^ 

1. As shown in |Sj, jJJ for special examples, f, in a larger sector is 
expected to have singularities with an accumulation point at infinity. 

2. In section IHl we also show that in some special cases, there is a duality 
between small t and large x. 

3. Relatively simple examples in which the assumptions apply after suit- 
able transformations are the modified Harry-Dym equation Hf + = 
H^Hxxx — Kuramoto-Sivashinsky ut + uux + Uxx + Uxxxx = and 
thin-film equation ht + 'S/ ■ (/I'^VA/i) = (the latter with initial conditions 
such as /i(x, 0) = 1 -|- (1 -|- axf + 6x2)"^ in d = 2). The former equation is 
discussed in detail in and the normalizing process, adapted to short time 
analysis, is described in ^ 

4. The condition a,- > 1 is not particularly restrictive in problems with 
algebraically decaying coefficients. For these, as discussed in ^8_, one can 



restriction of the form |x|°|r(x,t)| < Ar{T){*) may appear more natural. However, 
since every component of x is bounded below in 2?,(>,p(,,x, it is clear that (*) implies IIH with 
ftr ~ a/d. The same comment applies for condition CHJ- This form is more convenient 
in the present analysis. See also Note 4 following Theorem 
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redefine f by subtracting out from it the first few terms of its formal asymp- 
totic expansion for large x. The new f decays faster at oo and the condition 
to Or > 1 can be ensured. 

2.3. Borel summability of power series solutions and their asymp- 
totic character. Determining asymptotic properties of solutions of PDEs 
is substantially more difficult than the corresponding question for ODEs. 
Borel-Laplace techniques however provide a well suited modality to over- 
come this difficulty. The paper shows that formal series solutions are Borel 
summable to actual solutions (a fortiori are asymptotic to them). A few 
notes on Borel summability are found in ^7.21 

In addition to hypothesis of Theorem ^ we need, first of all, to impose 
restrictions to ensure that there exist series solutions, to which end the 
coefficients of the equation should be expandable for large x. In many 
practical applications these coefficients turn out to be finite combinations of 
ramified inverse powers of Xj. 

Condition 2. For large |x| and some N G N'^, the functions bq^k(x, i) and 
r(x, t) are analytic in {x-^^^^^ , ...,x^^^^'') 

Theorem 2. If Condition\^and the assumptions of Theorem^are satisfied, 
then the unique solution f found there is the Borel sum of its own asymptotic 
series. More precisely, f can be written as 

(14) f(x,t)=/ e~P-^^Fi(p,t)dp 

1 1 

where Fi is (a) analytic at zero in {p'^'^'^ , ...,p^ '*); (b) analytic in p ^ 
in the poly-sector |argpj| < ;^^^</' + 2{n-i) ' ^ — ^'^'^ exponentially 
bounded in the latter poly-sector. 

Comment: For PDEs it is known that it difficult to show, by classical 
methods, the existence of actual solutions given formal ones, when the formal 
solutions diverge. Borel summability of a formal asymptotic series solution 
shows in particular, using Watson's lemma [3], that there always indeed exist 
actual solutions of the PDE asymptotic to it. Borel summability also entails 
uniqueness of the actual solution if a sufficiently large sector of asymptoticity 
is prescribed (see, e.g., PP). The Borel summability parameters proven in 
the present paper are optimal, as explained in the following remarks, and 
the sharp Gevrey class of the formal solutions follows too. 

n 

Remark 3. (i) It follows from the same proof that x "-i can be replaced with 
x'^ for any (5 G [1, ;^^]. The canonical variable in Borel summation is that 
in which the generic Gevrey class of the formal series solution is one (i.e., 
the series diverge factorially, with factorial power one; ^J. This variable, 

n 

in our case, is x"-i . 

(a) At least in simple examples, the sector of summability is optimal. See 
also Note\4^ 
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(in) In many problems of interest the conditions of Theorem\^are met by 
the equation in more than one sector (after suitable rotation of coordinates). 
Then the functions Fi obtained in ^ are analytic continuations of each- 
other, as it follows from their construction. 

(iv) If we had made the change of variable x x^/i"'-'^) first, (yielding 
the normalized Borel variable), the transformed PDE would have been more 
difficult to handle. Borel transforming directly from the x to p instead re- 
quires us to perform, in the proof of Theorem\^ an acceleration in the sense 
of Ecalle to establish Borel summability, but is technically simpler. 

The proof of Theorem [21 is given in ^ 
See also H7.11 

2.4. Spontaneous formation of singularities in nonlinear PDEs. Borel 
summability of formal solutions associated to solutions with prescribed ini- 
tial data is a key ingredient in the detailed analysis of spontaneous singu- 
larities of solutions and in the study of their global properties. Applications 
of the present techniques in these directions, partly relying on extensions to 
PDEs of the methods in 0, are discussed in the paper [TT] . 



3. Inverse Laplace transform and associated integral equation 

The inverse Laplace transform (ILT) G(p, t) of a function g(x, t) analytic 
in X in T)^ p.y^ and vanishing algebraically as x ^ oo (cf. Lemma |1] below 
and Note following it) is given by: 

(15) G(p, t) = [£-i{g}] (P, t) = -i-j / eP-g(x, i)dx 

[2ttiY Jcd^ 

with a contour Cd as in Fig. 1 (modulo homotopies), C 'D^^p-:x:, and p 
restricted to the dual (polar) domain 5^ defined by 

(16) S^ = {p: \pi\ > 0; argp^ G {-(p,(p),i = l,...,d} 

to ensure convergence of the integral. 

The following lemma connects the p behavior of the ILT of functions of 
the type considered in this paper to their assumed behavior in x. 

Lemma 4. // g(x, t) is analytic for x in T>^^p.^, and satisfies 

(17) |x"| |g(x,t)| < A(r) 

for a > ao > 0; then for any 6 £ (0, (p), the ILT G = C^^g exists in S^ps 
and satisfies 



A{T) 

for some C = C(5, qq) 



(18) |G(p,t)|<C^|p"-i|e2|Pl'' 
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Proof. The proof is a higher dimensional version of that of Lemma 3.1 in 
[Hj. We first consider the case when 2 > a > uq. Let Cp^ be a contom' 
so that the integration path in each x component is as shown in Fig. 1: it 
passes through point pi + and s = pi + \pi\~^ + exp(i0signum(r)) 

with r £ (—00,00). Choosing 2p> pi> {2/y/?,)p, we have |s| > /3 along the 
contom' and therefore, with arg(pj) = 6 £ {—(p + d^cp — 6), 



;(s,t)| < ^(r)|s-"| and |e^-P| < c^iIpI+VIpI'^^^I'^+^'I 



Thus, 



e^'Pg(s,t)ds 



pi 



< 2^(r)e''i|Pl+'^JJ 



Pi + \Pi\ +ire 



-|pi|rsin5 



dr 



(19) 



< 



where K and K are constants independent of any parameter. Thus, the 
Lemma follows for 2 > a > ag, if we note that r(a) is bounded in this 
range of q, the bound only depending on ao- 

For a > 2, there exists an integer A; > so that a — k £ (1, 2]. Taking 



[{k - l)!]'^h(x, t) = r g(z, t)(x - z)'^-idz 
(clearly h is analytic in x, in T^^^p and (9]^h(x, t) = g(x, t)), we get 



^^"'') = I(f^A g(x.y,t)(y-l)(^-i)My 



( _-\\dk {k-a)l [-co 

^ P _!),]. / A(x . y, t)y-(y - l^^^-^^^dy 



with |A(x • p, t)| < A{T), whence 

|h(x t)\ < -^(^)[^("-^)]' 

' ^' - |xi|°-nr(a)]'^ 

From the arguments above with a — k playing the role of a, we get 



|£-Hh}(p,t)| < C(<5)^^^|pir-'=-ie2|p|P 



Since G(p,t) = (— l)'^'^p-'-'^£ ^{h.}{p,t), by multiplying the above equation 
by Ip"'"!'^, the Lemma follows for a > 2 as well. □ 

Remark 5. The constant 2p in the exponential bound can be lowered to 
p + 0, but m8\) suffices for our purposes. Note also that the statement also 
holds for p = 0, a fact that will be used in ^6. 
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Remark 6. Corollary\^ below implies that for any p G S^, the ILT exists 
for the functions r(x, t), bq^k(x, t), as well as for the solution f (x, t), whose 
existence is shown in the sequel. 

Remark 7. Conversely, i/G(p, t) is any integrable function satisfying the 
exponential bound in U8\) . it is clear that the Laplace Transform along a ray 
(PP exists and defines an analytic function of x in the half-plane for each 
component defined by ^[e^^'Xi] > 2p for Oi G {—(p,cl)). Due to the width of 
the sector it is easy to see, by Fubini, that CG = g. 



Remark 8. The next corollary finds bounds for Bq k = ^ ""^{bq^k} o.iT'd 
o/bq_k and r in V^^p^ D V^^p. 



R = £ ^{r} independent o/arg pi for p G Sa, following from the properties 



Corollary 9. The ILT of the coefficients bq^k (cf. ilO^) and of the inho- 
mogeneous term r(x, t) satisfy the following upper bounds for any p G 5^ 

(20) |Bq,k(p,t)| < ^^^A(r)|p"'>-i|e2''o|pl 



(21) |R(p,t)| < ^i^A,(r)|p--i|e2P0lPl 

[r(ar)]^ 

Proof. The proof is similar to that of Corollary 3.2 in 8 . From the conditions 
assmned we see that bq k is analytic in x G T)^^^p^yy^ for any 0i satisfying 
(2n)~^TT >(/>!>(/)> 0. So Lemma m can be applied, with g(x, t) = bqk, 
with 01 = + ((2n)^^7r — 0)/2 replacing (/>, and with 8 replaced by (/>i — = 
((2n)"^7r - (/))/2. The same applies to R(p,t), leading to and (ED). In 
the latter case, since > 1, ao in Lemma^can be chosen to be 1. Thus, 
one can choose C2 to be independent of a,.. □ 

Lemma 10. For some R G and all p with |p| > R and maxj<rf | argp,| < 
(j) we have for some C > 

(22) 5RP,(-p) >C|pr 

Proof. For the proof, we take S = {p : |p| = 1, maxj<(^ | argpj | < (j)} and 
note that 

(23) Co = inf m^n-A-V) > 

pe-B 

l<j<m 

(cf. definitions following ©). Indeed, if Co = 0, then by continuity 
^'Pn;ji—p) would have a root in B which is ruled out by @. The con- 
clusion now follows, since on a sphere of large radius R, Vj is given by 

R''Vn;j{-p/R) + 0{R^). □ 
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The formal inverse Laplace transform (Borel transform) of Q with respect 
to X (see also COI) for p E 5,^ is 



(24) dtF + r{-p)F = J^J2 * F*'' * ((-P)'^O + ^^P' 

where F = £~^f . After inverting the differential operator on the left side of 
(|24|) with respect to t, we obtain the integral equation 



(25) F(p,t) =AA(F) = Fo(p,t) 

+ re-^(-P)(*-^)5]'5]Bq,k(p,r)*F*'^(p,r)**n ((-p)^^^/ ( P, r))*'^''JdT 



q^O kbO l,\i 

where 
(26) 

Fo(p,t) = e-^(-P)*F/(p)+ / e-^(-P)(*-")R(p,r)dT and F/ = ^-^f/} 

JO 

Our strategy is to reduce the problem of existence and uniqueness of a 
solution of © to the problem of existence and uniqueness of a solution of 
((2SJ, under appropriate conditions. 

4. Solution to the associated integral equation 



To establish the existence and uniqueness in (|25j) we first introduce suit- 
able function spaces. 

Definition 11. Denoting by S,j, the closure ofS^ defined in dS^ = 

and /C = 5,^ x [0,T], we define for v > (later to be taken appropriately 
large) the norm \\ ■ \\u as 

(27) \\G\\. = M^ sup (ll{l + \pAe-'\P\\G{p,t)\ 

(p,t)G/C V i / 

where the constant Mq (about 3.76) is defined as 

, , f2(l + s2) (ln(l + 5^) + s arctans) 1 

28 Mo = sup ^ \ 2/., -\ 

s>o [ s{s^ + 4) J 

Note.- For fixed F, \\F\\i, is nonincreasing in u. 



Definition 12. Consider the following Banach space. 

(29) = {F : F{-,t) analytic in 

and continuous in for t G [0,T] s.t.||F||i, < ooj 
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Remark 13. // G G A,f,, then g(x, t) =: CelG} exists for suitable if 
pcos{6i + avgXi) > V. Furthermore, g(x,f) is analytic in x, and |x-'-g(x,i)| 
is hounded in V^^p-^. 

Lemma 14. For v > Apo + Ur, F/ in \2f^) satisfies 
while R satisfies the inequality 

\M\v < c(</')^r(r)(z^/2)-'^°'-+'^ 

and therefore 

(30) llFoll. < C{cl>)Ao{T){u/2)-'"'^+^ 

Proof. This proof is similar to that of Lemma 4.4 in |H]. We use H21|l . note 
that Oj. > 1 and also that for v > ApQ + we have 

(31) sup ^ll e-(-2po)|P.| < [ar±i)^.ar^i _ ^^^y^r^i 
|pi|>o r(ar) T{ar) 

where K is independent of v and a^- The latter inequality follows from 
Stirling's formula for T{ar) for large a^- 

Using the definition of the z/— norm and the two equations above, the in- 
equality for ||R||i/ follows. Since f/(x) is required to satisfy the same bounds 
as r(x,t), a similar inequality holds for ||F/||,y. Now, from the relation (|26() 
and the fact that 5R'Pj(— p) is, by Lemma^l bounded below for p G 5,^, we 
get the following inequality, implying (|3U() 



|Fo(p,t)| < |F7(p)| +r^o(T) sup |i?(p,t)| 

Q<t<T 



□ 



It is convenient to introduce a space of sectorially analytic functions pos- 
sibly unbounded at the origin but integrable. 

Definition 15. Let 

W := |h : H(p,t) analytic in 5<^,|H(p,t)| < C |p"-^| e''IPl} 

(C , a and p may depend on H). 

Lemma 16. // H G H and F G A^, then for u > p + 4, for any j, H* Fj G 
A,/,, anS: 

(32) ||H * FjW^ < |||H| * \Fj\\\^ < C[T{a)f 2'^"(zy - p)-'^°||F||^ 
where C is independent of a. 



^In the following equation, || ■ is extended naturally to functions which are only 
continuous in K.. 
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Proof. The proof is a vector adaptation of that of Lemma 4.6 in [Hj. From 
the elementary properties of convolution, it is clear that H * Fj is analytic 
in S(f) and continuous in 5^. Let 9i = aigpi- We have 



|H*F,(p)|<||H|*|i^,|(p)l< / |H(se*^)||F,(p 

Now 

(33) |H(se*^)| < Cls^-ijel^l" 

and 



(34) / s"-ie'"'^|Fj(p-se*^)|(is 

•''njo.ip.i] 



< ||Fj||,e^lP'|p-| J] 



Jo Mo(l 



'^-'^Q-i'^-p)\Pi\si 



+ IkP(i-s.)' 



-dsi 



Since v — p > 4, we can readily use p22jl in the Appendix with n = \pi\, v 
replaced by — /O, cr = 1 and m = 1 to conclude 

(35) ip.r /' .,f'' < - 



JO Mo(l + lBr^(l - Si)^) A/„(l + |pi|2) 

Therefore, from (|34p. we obtain 
(36) 

/ s-^el^hF,(p - se^«)|ds < ^[r(«)]^ h^.^!!' I "pf 

^njo,ip,i] ^^''^i(l + |p^r) 

From this relation, ()32() follows by applying the definition of || • □ 
Remark 17. Lemma [7^ holds for p = as well, when u > A. 

Corollary 18. For F G A,f,, and v > Apo + 4 we have Bq^k * Fi ^ -A-cj, and 

||Bq,k*Fn|. < |||Bq,k|*|F||L<TO(,/.,aq) (z./4)-'^°Mf,(r) ||F||, 

Proof. The proof follows simply by using Lemma 1161 with H replaced by 
Bq,k and using the relations in Corollary |3 □ 

Lemma 19. For F G Acj,, with v > ApQ + 4, for any j, I, 
Proof. From the definition ((JJ, it readily follows that 

|Bq,k*(pji^«)l < |pj||Bq,k| *|^/| 

The rest follows from Corollary ()18p. and the definition of || • □ 
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Lemma 20. For F, G G and j >0 

(37) \{piFi,)*Gi,\<\p^\ ||F|*|G|| 

Proof. Let p = (pie*^^,p2e*^^, ••,Pde*^'*). Then the result follows from the 
inequality 



(38) 



s>Fi^{s)Gi2{p - s)ds 



< |pJ| / |F(e^^s)||G(p-e*^s 

njo.ipi.] 

□ 



Corollary 21. IfFe As, then 



(39) 



< 



Proof. This follows simply from repeated application of Lemma I2UI □ 



Lemma 22. For F, G G As, 



IFI * IGI 



< 



|F|L||G|L 



Proof. 
(40) 



|F| * |G| 




/ 






Jo 



|F(s)||G(p-s)|ds 



< / |F(e^^s)||G(p-e^^s)|ds 

•''njo.ipi,! 



Using the definition of || • W^^, the above expression is bounded by 

dsi , Ipile'-IPl 



iiFiuiiG|i. [[j^ jrT:sf)[i + (Iki - s.r] ^ M^u^a + \p^\') 



< 



IFIIJIG 



The last inequality follows from the definition (|28|) of Mq since 



\pi\ 



dsi 



{l+sUl + {\p^\-sin 



,ln(|p,|^ + l) + |pi|tan-^|pi| 

\p^\{\V^? + ^) 



□ 



Corollary 23. For F, G G A^, then 

|||F|*|G|||^< ||F||^||G||^ 
Proof. This is an application of Lemma EH and the definition of 



□ 
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Lemma 24. For u > Apo + 4, 

(41) 



< g ILIP^I ||F||N+l'^l-i |Bckl*|F| 



if (q, k) 7^ (0, 0) and is zero if (q, k) = (0, 0). 

Proof. For (q, k) = (0, 0) we have Bq^k = (see remarks after eq. (|10|) ). If 
k 7^ 0, Corollarv 1211 shows that the left hand side of 1)411) is bounded by 



1*9;, 



|Bq,k| * |F| * |F|*(l''l-i) * Y{\F 

Using Corollaries |^ and ESI ^-nd Lemma 1221 the proof follows for k 7^ 0. 
Similar steps work for the case k = and q 7^ 0, except that Bq k is 
convolved with pj -Fi^ for some (j', h), for which the corresponding g/^ j/ 7^ 0, 
and we now use Lemma QUI and the definition of || • □ 



Corollary 25. For v > ApQ + 4, 



(42) 



Bq,k*F*'^* 11(1^^' 



< 



|kl + |k| 



The proof follows immediately from Corollary and Lemma □. 
Lemma 26. For v > 4po + 4, we have 



dr 



^ cA(r)e-ipi^ rr'"^iFiii''i+i-i 



(43) 



for some Af,(T) > Af,(T) (evaluated in the proof) and where the constant C 
is independent ofT, but depends on (p and a^. 

Proof. This is a consequence of Lemmas and [51] and the fact that for 
< |1'| < n we have, for |p| < R (with R as in Lemma ITU]) . 



(44) 



J := 



P 



Jo 
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For IpI > Rwe have, by Lemma ITTH V{—p) > C|p|", and J is majorized by 



(45) mmax ^ 



1 _ g-3?T',(-p)t 



Ti-|i'l/"|p|l''l 1-e- 



■7 



^ max — — sup 



j<m RVj{-p) L J - j<m |SRp^.(_p)||l'l/n 

where 1' = Xlj.J^ij- □ 

Definition 27. For F and h in A^, and Bq ^ G "H, as above, define ho = 
and for k > 1, 

(46) hk = Bq,k*[(F + h)*'^-F*'^]. 

Lemma 28. For v > 4po + 4, and for k 7^ 0, 

(47) llhkll. < |k|(||F||, + ||h||,)''''~'|||Bq,k| *|h|L 
and is zero for k = 0. 

Proof. The cases |k| = 0, 1 fohow from the definition of ho and (|46l) respec- 
tively. Assume formula (|47|) holds for all |k| < I. Then all multiindices of 
length / + 1 can be expressed as k + e^, where e^ G M"* is the m dimensional 
unit vector in the i-th direction, and |k| = /. 

||hk+e, II. = ||Bq,k * {Fi + hi) * (F + h)*!^ - Bq,k * Fi * F*||, 

= ||Bq,k*/ii*(F + h)*l'+Fi*hk||. 

Using (|^ for |k| = Z, we get 
< |||Bq,k| * |h|||j, (||F||i, + ||h||,,)' + /||F|| 

!^ fll-^lll^ 

+ ||h||^) |||Bq,k| *|h|L 

< (Z + 1)(||F||, + ||h||,y|||Bq,k| *|h|L 
Thus (jUI) holds for |k| = / + 1. □ 

Definition 29. For F G and h G A(f), and Bq.k as above define go = 0, 
and for |q| > 1, 



(48) gq = Bq,k* n(p'[^'+^']) ''-Bq,k* 11(1^^' 



Lemma 30. For v > 4po + 4, go = and for |q| > 1 



(49) |gq| 



< 



p 



e^'IPllql / \|q|-i 

"|F||,. + ||h||,. |||Bq,k|*|h| 



Mo'n.(l + lKp) 



I V 



and is zero for q = 0. 
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Proof. The cases |q| = 0, 1 follow from the definition of go and (|48j) respec- 
tively (since only terms linear in F are involved in (|48j) ). Assuming (|49|) 
holds if |q| < I we show that it holds for q + e, where e is a unit vector, say 
in the {h, j'l, j'2j j'd) direction. We have 



(50) |gq+e|< 



Bq,k * 



^'{Fi,+hi,)\ * *ll [p^{Fi + hi) 



Bq,k * 



< 



Bq,k * (P' /l/i 



pj>/i] * *n [p"^^ 



*qi, 



+ l(P''^«l)*gq| 



Using Lemma OH and equation ()49() . we get the following upper bound im- 
plying the induction step 



F||^ + ||h||^)^'?''J |||Bq,k|*|h| 



+ 



< 



|q|e' 



|pEjte,j+ei,j)(|q + e|e^|P 



F||, + ||h||,)l'il-i||F||,|||Bq,k| *|h||L 
F||^ + ||h||^)l'il|||Bq,k| *|h|L 



□ 



Lemma 31. For F and h in A,j,, v > 4po + 4, 

*k 



Bq,k * (f + h) '"'^ * (pj {Fi + /i;)) - Bq,k * F*!^ * II {j^Fl 



IpEj^ijI (|q| + |k|)e^lPl 



^^^^ " Mo^n.(i + iKP) 

if (q, k) 7^ (0, 0) and is zero otherwise. 



F||, + ||h||,)l'^l+l'^|-i|||Bq,k|*|h| 



Proof. It is clear from ()46() that the left side of (|51() is simply 
hk* *n (p'(i^i+/i/))*'''+F*'^*gq 
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However, from Corollary I21L Lemmas 1221 and I28 ( 



and from Corollary 1211 Lemmas 1221 and I3U1 



F||, + ||h||,)l'^l + l^|-M||Bq,k| 



< 



P 



|q|e 



''IpI 



|F||, + ||h|L)l''l+l'il-Ml|Bq,k|*|h 



Combining these two inequalities, the proof of the lemma follows. 



□ 



Lemma 32. For v > 4po + 4: we have 



(52) 



-P(-p){t-r) 



*k 



dT 



< A(r)c(</.)(|q| + |k|) m\u + ||h|U)W+i^i-i 



-daa 



|h| 



Proof. This follows from Corollary 1181 and Lemma 1^ and the definition of 
II • III, together with the bounds and (|^ . □ 

Lemma 33. For F G A^, and v > 4po + Or + 3 large enough, (see Note 
after Definition ), J\f{F) defined in \2S\) satisfies the following hounds 

(53) ||AA(F)|U < llFoll. + C(<A)i.(r)^' Q'''^" ||F||l'il+W 



(54) ||Ar(F + h) - AA(F)||,. < C(,^)i,(r)||h||,x 



-da„ 



(|q| + |k|)(||F||,.+ ||h|U)l'il+W-i 



Proof. The proofs are immediate from the expression 1)25^ of A^(F) and 
Lemmas 1261 1281 and 1321 Note also that the sum with respect to q only 
inyolyes finitely many terms, see ©• 

□ 



Remark 34. Lemma ^SM is the key to showing the existence and uniqueness 
of a solution in to \25^ . since it provides the conditions for the nonlin- 
ear operator M to map a ball into itself as well the necessary contractivity 
condition. 
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Lemma 35. // there exists some b > 1 so that 

(55) 6||Fo||. < 1 
and 

(56) C7(c^)A(T)E'E WbFot^^^-^ <l-l 

then the nonlinear mapping J\f , as defined in \25)) . maps a hall of radius 
6||Fo||i/ into itself. Furthermore, if 

(57) C(<A)i6(r)j;' J^dql +|k|) (^)"'"''(36)l'^l+l'il-i||Fo||l'^l+l^l-i < 1 
then M is a contraction there. 

Proof. This is a simple application of Lemma lHHl if we note that in the ball of 
radius &||Fo||, ||F||^ < 6'^||Fo||^ and using in ((531) t^i^ fact that ||F||,^ + ||h||;^ < 
3611 FolU if inax{||F||j,, ||F + h||i^} < 6||Fo||. □ 



Lemma 36. Consider T > and G (0, {2n)^^'ir) so that (0) is satisfied. 
Then, for all sufficiently large v, there exists a unique F E Acj, that satisfies 
the integral equation Ii2,5\). 

Proof. We choose 6 = 2 for definiteness. It is clear from the bounds on 
||Fo||i/ in Lemma El that for given T, since a,- > 1, we have 6||Fo||i/ < 1 for 
all large. Further, it is clear by inspection that all conditions (|55|) . (|56|) 
and (|57j) are satisfied for all sufficiently large ly. The lemma now follows 
from the contractive mapping theorem. □ 

4.1. Behavior of T near p = 0. In the following proposition, we denote 
by T the solution F of Lemma l36l 

Proposition 37. For some Ki > and small p we have \^\ < i^i|p-'-|"''~^ 
and thus |T| < i^2|x"'"|^"'' for some K2 > in T^^^p as |x| — > 00. 

Proof. The idea of the proof is to note that, once we have found T, this 
function also satisfies in a neighborhood of the origin Sa = ^njp : < fli} 
a linear equation of the form 

(58) T = g(T) + Fo or T = (1 - g)-iFo 

where, of course, Q depends on the previously found T; there are many 
choices of Q that work. Every term in the sum in (|25|) is a convolution 
product; in each of them we replace all but one component of F by the 
corresponding component of ®F; QY is defined as the sum of the terms 
thus constructed. Estimates of the form used for Lemma 1551 show uniform 
convergence of the sum for large enough v (or small a). The result is a ^ as 
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below, where the sum over fi contains only finitely many terms and which 
has manifestly small norm if a is small (or i' is large) 



(59) 



GF 



-P(-p)(t-r) 



Y^Gi*Fi + Y,G^*{{-prFi^] 



dr 

I in 5« 



By (dU, ((EI), ^ and LemmalU we see that ||Fo||oo < K3 |a° 
for some > independent of a. Then, from (|58|1 for small enough |a|, 
we have 



m_ax|'F(p,t)| = ||T|| < (1 - ||g||)^^max||Fo|| < 2Ks\a.''--^\ 

Sa Sa 

and thus for small |p|, we have |F(p,t)| < 2X3 |p°''^"'^| and the proposition 
follows. Indeed, the arguments also show that that the same estimates hold 
when any component pi 0, if the others are bounded. □ 

4.2. End of proof of Theorem^ Lemma^ shows that if f is a solution 
of © satisfying |xi||f| < A{T) for x G 2?^,p,x, then £~i{f} G A^s for 
< 6 < (J) for v sufficiently large. For large enough p, the series (|lfl|) 
converges uniformly for x G ^</>,p,x and thus F = £~^{f} satisfies (|25|) . 
which by Lemma l36l has a unique solution in for any (p G (0, (2n)"^7r) 
for which Q holds. Conversely, if T G A^ is the solution of ()25() for 
u > ui, then, for sufficiently large = £ T is analytic in x in I'<^,p 

for < (f) < (j) < (2n)~^7r (cf. Remark [T^ . Proposition IS7I shows that 
^ = 0(x~"'') and entails uniform convergence of the series in (jSJ. By the 
properties of Laplace transforms, T solves the problem (jSJ. 



5. BOREL SUMMABILITY OF FORMAL SOLUTIONS TO THE PDE 

We now assume Condition 1 in addition to Assumption 1. In our ap- 
proach it was technically convenient to use over summation, in that the in- 
verse Laplace transform was performed with respect to x. Showing Borel 

n 

summability in the appropriate variable (x"-i, as explained) requires fur- 
ther arguments. 

5.1. Behavior of F for large |p| outside 5^. For the purpose of showing 
Borel summability of formal series solutions we need to control F for large |p| 
uniformly in C^. For this purpose we introduce two other Banach spaces, 
relevant to the properties we are aiming to show. Firstly, let !B(z^, n,5) 
be the Banach space of functions analytic in the sector 5 = {p : \pi\ > 
0, arg(pj) G {ai,bi)} and continuous in its closure, where bi — Oj will be 
chosen larger than 27rNi (cf. Condition ^ The Banach space is equipped 
with the norm 

(60) \\^\U= sup *(p,t)e-^(*+i)S,(fel+IP.r) 

peS;telo,T] 
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Lemma 38. For any intervals {ai,bi), i = l,...,d the solution F of \25\) 
given in Lemni,a [9^1 is in Q3(z^, n,5). 

Proof. Because of the obvious embeddings, it suffices to sliow that for any 
S, (|25|) has a unique solution in *B(z/, The proof of this property is 
very close to that of Lemma l36l after adaptations of the inequalities to the 
new norms, which are explained in the Appendix, ^7.41 □ 

5.2. Ramification of F at p = and global properties. We define 
*B(z^, n,ei) to be the Banach space of functions defined on Sf_^ = {p : 
maxj \pi\ < ei} in the norm (Ififl with S replaced by Sf^. 

Lemma 39. Let 

h 2d_ 

(61) G(p)= Pi'-'-Pd'An^.^Ap) 

where ^^1,...,^^ are analytic at p = 0. Then the functions are unique 

and for some constants Ci and C2 and large p we have 

(62) \Aj,_j,ip)\ < Ci\pf' max \G{pie^^^^\ ...,Pde^^''n\ 
In particular, in Sf we have, for some constants C3 and C4, 

(63) C3 max sup \G{pie'^^'''\ ...,pde^^''''')\ 

< sup |^,„...,,,(p)| < C4 max sup \G{pie^^'^\ ...,Pde^^'^^')\ 

Remark 40. We note that in \6'/^) the order of analytic continuations is 
immaterial. 

Proof. The proof is by induction on d. We take d>\, assume (|39() with A-^ 
analytic and write p = {pi,p^). We have 

(64) G{p)= Yl Pi'i E P2'---P^'An,-,uiP) 

0<ji<Ni \{i™.<Af™;m=2,...,d} 

=: Y pf'G.M.P^) 

0<ji<Ni 

(with the convention that Gj-^ = Aj^ if d = 1). We write the system 
(65) 

G{pie'^^\p^) = Y e2'=WA^Vi^G,^(pi,p^); A: = 0, 1, iVi - 1 

0<ji<Ni 

which has nonzero Vandermonde determinant, from which Gj^{pi,p^) are 
uniquely determined, which in turn, by the induction hypothesis determine 
^iivjd' with the required estimates. □ 
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Lemma 41. Under the assumption 1 and condition 1, the solution in Lemma 
1. 961 can he decomposed as follows: 

(66) F(p,t)= pf^ ■■■pf'HP,t) 

0^j-<N 

where Aj(p,t) G !B(z^, n,5) are analytic atp = 0. Furthermore, in analyzing 
the continuations in restricted sectors pe^'^*^ G 5^ we have for some u, in 
the norm defined in \2T\ ) (cf. also Remark \4(A ) 

(67) max{||F(-e2-^J,.)|U,{||Aj(-,-)||.};0^j -<n} < OO 

Proof. We consider the equation (|25)) on 53(z^, n,5)^ where N counts the 
Aj(-,t) via the decomposition (|66|) . Noting that 

^ ^ ^ r(a + /3 + 2) ^ 

it is straightforward to show that the space of functions of the form (|61() 
is stable under convolution. Since R{p,t) and therefore Fo(p, t) are of the 
form ()66() it follows that M leaves the space of F of the form H66|) invariant. 
Using the estimates (|63() we see that M is well defined in a small ball of 
radius £2 in in *B(i/, n,5) and that it is a contraction there. Therefore the 
solution to (Uni) is of the form For pe^'^^J e 5<^, ||F(pe2'^*j)||,, are 

well defined. Using again Lemma OH the first statement follows. To show 
finiteness of || Aj(-, t)||i, it suffices to prove finiteness of ||F(pe2''*j)||,,. To this 
end, we note that all these continuations satisfy equations of the type © 
with coefficients satisfying the requirements in ^and thus the result follows 
from Lemma 

□ 

Lemma 42. Assume G is an entire function of exponential order n, more 
precisely satisfying the inequality |G(p)| < Ce'^'p'" for some constants C,u 
and that in a sector 5^ = {p : |p| > 0, maxj | arg(pj)| < (j)}, it grows at most 
exponentially, |G(p)| < Ce'^^'P'. Then there exists a function Gi increasing 
at most exponentially |Gi(p)| < Ce'^^'P' in any proper subsector of S^-^ where 
01 = 2{^n-i) ~'~ ^^'^ such that G(z") is analytic at z = 0, such that 

(69) g(x):=/ e-P-^G{p)dp= e^P'^-^ Gi(p)dp 

Jo Jo 

Proof. We start with the case when G, x and p are scalar, the general case 

following in a quite straightforward way as outlined at the end. 

The assumptions on G ensure that the first integral in (|69|1 exists and 

g{x) has an asymptotic power series in powers of x~^ in a sector of opening 

7r + 2(p centered on R"*". The function gi{x) = g{x^'^~^^^^) has a (noninteger) 

power series asymptotics in a sector of opening ^^(tt + 2(j)) and by the 

general theory of Laplace transforms, Gi := C^^gi is analytic in a sector 



NONLINEAR EVOLUTION PDES IN R+ x 23 

of opening (vr + 20) — tt centered on M^, Laplace transformable, with. 

Laplace transform gi. It follows that 

(70) 

1 pc+ioo poo poo 

Glip) = ^ e-'^^' "''G{q)dqdu =: / K^{p,q)G{q)dq 

Jc-ioo Jo Jo 

We show that Gi has a convergent expansion in powers of p^/" at zero. The 
function 

(71) {p, q) = (^) " {qVp^-') 

is Ecalle's acceleration kernel ^ E|- For a G (0,1), with /3 = 1 — a, 
c = iSa"^^, the function Ca is an entire function and has the following 
asymptotic behavior pi I14j: 

1 

(72) Ca{x) ~ -^==xi/2g-cx. ^ I g^^g^i ^ I 
Using ((7T|) we see that 



(73) 



roo poo 

/ K^[p,q)q^dq=p^''^~^^^^/'' s'^+'^Cn^ (s")ds 



We expand the entire function G in series about the origin, G{q) = XlfeLi^ 9kQ^+ 

Rn^q) and note that 

(74) 

oo oo 

\RN{q)\ < Yl IG^'^^WIkl'A! < IC^'^^WIkl'A! < C'e-^'"'" = E{q) 
k=N k=0 

uniformly in C. By JZ^ and ^ E{q)Ga{q'' /p^'^) is, for small enough p, 
in Li[0, oo] in q. By dominated convergence, we have 

j-oo poo 

/ K n^{p, q)G{q)dq = lim / Kn^{p,q) Y gkq'^dq 

Jo " N^OoJq n ^ 

and, using (|7H|) it follows that for small p, Gi is the sum of a convergent 
series in powers of p^/", as stated^. 

The argument for d variables and vectorial G is nearly the same: a vec- 
torial G is treated componentwise, while the assumptions ensure that the 
multidimensional integrals involved can be taken iteratively, the estimates 
being preserved in the process. 

□ 

Collecting the results of Lemma ^2 and Lemma applied to each of the 
Aj, the proof of Theorem [21 follows. 

■^To estimate the radius of convergence of this series it is convenient to start from the 
duahty (1691 and apply Watson's lemma, using Cauchy's formula on a circle of radius 
k^^"/{nuy^" to bound |G('='(0)| . 
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Note 43. In the example dtu+ {—dx)^u = we have = Formal expo- 
nential solutions have the behavior, to leading order, exp (^Cn{—x) "-i t^n^ 

n n 

with Cn = {n — l)/4/n""i (for all determinations of (— ). This also 

71 

points to as natural variable and indicates that the sector of summa- 

bility cannot be improved since it is bordered by (anti) stokes lines. 

6. Short time existence and asymptotics in special cases 

In some cases, the Borel summation approach can be adapted to study 
short time existence of sectorial solutions and study small time asymptotics. 
One important application is in the analysis of singularity formation in PDEs 
For simplicity, and since some assumptions are less general than in the 
rest of the paper, we restrict to d = 1 (scalar case) in this section. 

We motivate the assumptions made by looking at a particular example 
arising in Hele-Shaw flow with surface tension 

(75) Ht = -^ + H^H,,,, H{z,0)=z-'/^ 

the modified Harry-Dym equation (see > where it arises with ^ = z + 1 (as 
a local approximation near an initial zero of the derivative of a conformal 
mapping) . 

6.1. Formal series, preparation of normal form. Note: To simplify 
notation, in the following we let p stand for generic polynomials, for 
polynomials with nonnegative coefficients, and p(„) for polynomials of degree 
n. Similar conventions are followed for \} which represents homogeneous 
polynomials. Substituting in (|75|l a power-series of the form X^J^Lo ^"-^"•(■^) 
where Hq = z"^/^ yields the recurrence 

(76) ^ 

{n+l)Hn = -- ^ Hn^Hn2Hn^+ ^ Hn^Hn2Hn^H'^^ 

which inductively shows that Hn = z^'^^'^\]^^'^{z^'^/'^ , z^^). We let 
(77) 

N N 

gN{x,t) ■.= Y,t'"Hn{z) = J^f)(„)(te-3^te-2/3); where x = -z^'^ 

k=0 n=0 

In terms of x, (|75|) becomes, 

1 3 1 

(78) ATiH) := Ht + -H^^ - y^^H^-- " ^H^H^x + = 

It is straightforwardly shown that 

(79) MgN{x, t) = r ix-5p(4jv+i) (^x-^ tx'^/^) 

where for small xi,X2 we have moreover 

(80) P {AN +l){.Xl,X2) = \){N+1){X1,X2) [1 + 0(xi,X2)] 
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It is then natural to substitute : 

(81) H{z{x),t)=gN{x,t) + x-'^f{x,t) 

into (|75|) : we choose without loss of generality > 3. 

It will follow from the analysis that \ f{x,t)\ = o (j;^/^f)(Ar) (tx"'^, 
for smah t^/'^x^^ with argx G (-f - </>, f + </>) and (j) £ (0, |), thus ~ 
Yln=Qt'^Hn{z) for small t^/^x-i (see Corollary El) . 

Substitution shows that f{x,t) satisfies an equation of the form with 
n = 3 (third order, m = 1 (scalar case), with (cf. also ()1U() . and l|114|) below) 

(82) 

r(x,i) = t~ix^/=^P(4^+i)(te-3,te-2/3); b^f^ = x-^^'^x-"'i>*=Pq,fc;j(tx-3,te-2/3) 

Note: By (|8()|) . r(x, t) is small for small t or large x, in spite of the prefactor 

t-^x^/^ 

6.2. More general setting. Setting 1. We take po = 0, suitable for 
algebraic initial conditions in the domain, and consider the domain "D^fi^x, 
with 4> < small enough to ensure Q. Taking f(x,t) — f/(x) as the 
unknown function we may assume 

fiix) = 

(see Note 3 after Theorem and require that 

(83) |r(x,t)i <t-i^|xr^[)+,) (t^M^r'^s-.i^^'kr^^) 

where the degrees n'j satisfy 

(84) n'jPi - ujj > 1, for 1 < / < ET, 1 < j < Jr 

(As before, ()84() implies that r(x, t) is small for large x or small t). The pos- 
itive constants uJi, uj2-, ■■■^^Jri Pii P2, ■■■,Pk and 71, 72, ■■■,^k, are restricted 
by the condition 

(85) h := — > n 

71 

The labeling is chosen so that 

(86) n = — > — .... > — 

71 72 7i^ 

Also, if for some 1 < j < K — 1, ^ = we arrange f3j > The ujj 

are arranged increasingly: 

(87) oji < UJ2 < .... < wj^ 
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Furthermore, for any x E "D^^q^x, we require 

(88) |bq,k(x,t)| < {t^'\x\-^\.:;t^''\x\-^'') 

(89) /? > , aq,i > aq,2 > ... > aq,jq ; bq,k / ^ aqj + /3|k| > 
If only finitely many bq^k are nonzero we allow 

(90) /? > 

We also require that for all q, k for which bq^k 7^ we have 

(91) m-q^k := n + u;i(|q| - 1) - aq,i + {uji - /3)|k| - - ^jqij > 

Note: Assumption (|9H) is satisfied by modified Harry-Dym and by cer- 
tain classes of nonlinear PDEs and initial conditions- for instance, the thin- 
film equation ht + {h^hxxx)x = 0, with singular initial condition h{x, 0) = 
x~" for a > 0, but is generally quite restrictive. Weakening it requires more 
substantial modifications of the framework and will not be discussed here. 

Setting 2. Better properties are obtained under the assumptions de- 
scribed below. 

h = n 

V{-s) = 

(92) 

r(x,t) = \ Y.f=ix'^'^3 (tTij;-^i,....,t^^x-^^) 

bq,k(2;,t) = x-^l'^l Y^ili x-""^'^ a^,\^,j (rix-^i, ....,tT^^x-^^) 

where a^, ciq^kj are analytic near the origin and for small |z| we require, with 
the same restriction (|84|) on n^-, 

(93) |a,-(z)| <{)+,)(|zi|,...,|z„|) 

The restrictions on the numbers /32, ■■■Pk-, 71; 72, ■■■Ik-, Oqjj etc. are as 
in Setting 1. Furthermore, we assume that there is an u G so that the 
nonnegative numbers 
(94) 

"iq,k, ^^2-^1, -wi, aq,i-aq,2, ...,aq,i-aq,jq, ^72-/32, ..•,ra7x-/3x 

are integer multiples of nco. This condition, satisfied for the problem H75|) . 
comes out naturally in a number of examples and ensures the existence of 
a ramified variable in which the solutions are analytic. We choose > to 
be the largest with the property above. Define 

(95) C = yt~'/" ,f(c,t) = f(t^/"C,t) 
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and 

(96) D^,p = {C : \C\ > p; |argC| <<A} 

Theorem 3. (i) In Setting 1, under Assumption 1, there exists for large 
enough p a unique solution i{xt~^/^ ,t) to |3J), for ( = xt~^^''^ G '^(t>,p ^^'^^ 
with n'j as in i84\ ), 
(97) 

if(c,t)i < Elcr^^"^/"^l(n;)(lcr^^^^^-^^/"lcr^^•■•,^^"-^"/^cr^^ 

(a) In Setting 2, under Assumption 1, for any T > there is a p = p{T) > 
so that the mapping 

is analytic in V,),^p x {6 ■.\9\ <T}. 

Notes: 1. The function p will, generally, increase with T. 

2. The restriction (i = 1 is not essential, but made for the sake of simplic- 
ity. 

3. In these settings, there is a duality between large x and small t in the 
asymptotics: C, can be large either due to largeness of x or smallness of t. 
For f in a fixed interval, there exists some p so that the asymptotic bounds 
are satisfied for C, G D^ p. 

4. The following example shows that the requirement n > n is natural. In 
the equation gt + i—dx)^g = with g{x, 0) = substituting the expansion 
g{x,t) = X"" + Z]neN*"5n(a;), we get gn{x) = 0(x""""'). Thus one of the 
scales that emerge in the formal expansion is t/x"". On the other hand, in 
view of H83|) and (|88|) the most singular term as x — > is of the order t/x" 
since h = —. Combining with the above discussion we see that h > n. 

5. The leading order term in the Taylor expansion of 9 "^f, fo, satisfies 
an easily obtained ODE. The convergence of the series in part (ii) implies 
that singularities of fo can be related to actual singularities of the PDE for 
small time and this is the subject of another paper 

Corollary 44. For the initial value problem IjVS^ , for any T > there is a 
p = p{T) such that 

oo 

(98) H{z, t) = Y, t^Gkizt-"^'^) 

where the series converges in the region {{z,t) : |t| < T,\z\ > p, |arg2;| < 
Ivr} and Gk{C) are analytic in the sector {( : \(\ > p, \ arg("| < |vr}. 

6.3. Proof of Theorem 13] (i). It is convenient to make rescalings of vari- 
ables in Borel space as well. We note that 

/•oo 

(99) f((^,t) = ri/" / e-"^F{s,l;t)ds 

Jo 
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where 

(100) s = pt^/^ , F(s, A; t) = F(t-i/"s, tX) 

We use similar rescaling to define R(s, A; t), Bq^k(s, A; t) and Fo(s, A; t) where 
now 

(101) Fois,X;t) =tX e-*^^(-^*"'^")(i-")R(s,Ar;t)(ir 

Jo 



We let /iq,k = 1 - (^|q| + |k| + J2]=i TT=ij%j)- Using straight- 
forward calculations show that 

(102) ¥{s, A; t) = M{¥){s, A; t) = Fo(s, A; t) + J^' At^'^-x 



^0 



X , g -iAP(-st-i/^i)(i_^) 



m n 



1=1 i=i 



With slight abuse of notation we drop the hats from the newly defined 
functions. Let now 

(103) S^ = : args G {-(pA), < |s| < oo, < < ^| 

and consider the Banach space A(j, of analytic functions in 5^, continuous 
in Sfjj in the norm 

(104) ||F(-,-;t)||. = sup (l + |s|2)e-^l^l|F(s,A;t)| 

0<A<1, sG5^ 

Lemma 45. With r{x,t) satisfying i8!^) we have 

Jr 

||Fo(.,.;t)|U < e"*^Z^'^^+4('^^+^)/"{)+^. (^zy-^l,t^2-/32/nj^-/32^_^_^^7if-fe/nj,-/3^ 
J=l 

for V large (independent oft for small t), where —a is the lower bound of 

Proof. ^From ()83p . ()84|) and applying Lemma (with p = 0; see Remark 
we have 

|R(s,A;t)| 

tA "j V / 
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For A G (0, 1) we have e-*'^^-"*"'^")^^^-^) < e''* and thus (cf. THm ) 
(105) |Fo(s,A;t)| 

Jr 

i=i ^ 
Bounding each term of the polynomial in 11 • |L we obtain 

The proof now follows, choosing v sufficiently large and using (|84() and 1)86^ . 

(jHZI). □ 

Lemma 46. For large u, we have 

||Bq,k * Fill, < Cq,k(i^,t)||F II,,, where 

(106) 

co,o = 0; Cq,k(i^, t) = y-mti^-m)/n ^ K^.j^-^q,. t-«q.. /« ((q^ k) / 0) 
wzt/i constants independent of q, k, and t. 

Proof. Note first that bo,o = hence co,o = 0- From H88|l and Lemma 
(with p = 0), 

|Bq,k(p,t)| < IPl^i^i-'El^'l^'^Pk.- (^^Mp|^^^''^H^^•••,^^"|p|^") 

Switching from (p, t) to (s. A; t), 
|Bq,k(s,A;t)| <t(i-'^W)/"|s|^W-i 

X^|s|"<i.^t-°'i^^/>+^^. (^A^i|s|^i,A^2i72-/32/n|s|/32^___^^7A-i7K-fe/n|s|fe^ 

i=i 

For large z/, using Lemma (with p = 0) to bound in norm the terms of 

(107) ||Bq,k * F| < ||F||^ i(l-/3|k|)/n|^|-/3|k| 

J, 

|l/|-°'i.^t-"<i'^/"p+^^. f A^V-/^!, A^2^^2-/32/n^-/32^ ;^7Ki7A-fe/n^-fe 

Clearly, for large v, p^j^ can be replaced in (jlOTp by a constant Kj. Using 
l|86|) and (|89j) the conclusion follows. □ 
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Let now 



C{(f>, T) = max ■ 



I in \ I _ g-7 

peS^,\p\>Rfi<i'<n,^>o \ [-P) / 7 ' P&s^,\p\<Rfi<i'<n 
where R is the same as in the proof of Lemma [Till 

Lemma 47. For v large enough, M is contractive, and thus there exists 
unique solution F of ilOS^) . 

Proof. For v large enough, ()9ip . Lemma and Lemma HBl imply 
(108) C(</),r)^'^t'^-''Cq,k(i^,t)||2Fo||l'^l+l''l < llFolU 



and 



(109) C(</>,r)j;'j;t^-'<Cq,k(z^,t)(|q| + |k|)||6Fo||l'^l+l^|-i<l 

qbO k^O 

Now, Lemma (with /jq = 0, d = 1 and s replacing p), and Lemma EH 
imply 



m n 
— r *T 



- Mo(l + |s|2)'^'i''^^'^'*^ 11*11'^ 



Also, note that if /' > 0, s G with |st~i/"| > 
(110) 



< A 



1 - e 



The definition of C((/),r) implies that for I' > 0, s e with < R 

we have 



(111) 



^r^^-tP(-.t-V.)A(l-.)^^ 



Setting r = '^jqij, using (|110|) and (|111|) . we find after time integration 
(112) 



Ae 



-t-p(-st"l/")A(l-r) 



m n 

*-i — r * 



Bq,k*F*'^* n n(^'^0*''M^,AT;t)dT|U 



<//"-''/"C(</.,r)cq,k(^^,t)||F 



|q|+|k| 
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Using (jSU), (mU), (|TnH|) and (fTT^ , it follows that Af maps a ball of radius 
2 II Folio into itself. Using Lemma ISTl (|110)) and we obtain 



/ ABq,k*<^(F + h)** n Yl{s'[Fi. + hi]y 

I i=l 7 = 1 



m n 



1=1 i=i J 
< t''/"-'>C(0,r)(|q| + |k|)cq,k(i.,i) (||h||, + ||F||,)l'il+l'^l-i ||h|U 
where /' = Y^jqij from which the conclusion using H1U6|) and 1)91^ . □ 

Behavior of T near s = 

In the following proposition, we denote by T the solution F of Lemma 

113 

Proposition 48. For small s we have 

Jr 

1^1 < V \sr'^''h^^+'^^'^/^i)-^, (|s|^Si^2-/32/n|g|/32^ _^_^7k-/3xM|^|/3k) 

i=i 

Proof. The proof is similar to that of Proposition EZI using (|in5|l , (j83|) and 
(jHH). T to (|102j) solves a hnear equation 

(113) T = g(T) + Fo or T = (l-g)-iFo 

with ^ very similar to that given in §4. □ 

End of proof of Theorem 13] (i) The proof is a direct application of Lemma 
|17|and Proposition |1HJ Using and properties of Laplace transform, 
follows for large |C|, in the sector arg (— | — (/>, ^ + (/)). 

6.4. Proof of Theorem 13] (ii). An important difference is that infinite 
sums appear in some estimates. Analyticity of the functions a and the 
estimate 



pa-1 



-a+1 



Via) 

for V > 1 with C is independent of a and show convergence of the cor- 
responding series. Also, the proof of Lemma Wl\ holds if the following norm 
was used instead: 

||F||^= sup (l + |s|2)e-'^'^l|F(s,A;t)| 

0<A<l,|t|<T,sG5^ 

since for n = n, ^tV{—st~^^"') = is independent of t in the exponent 
in (Unp. To show analyticity, we let G{s,X;e) = 6-^^+"^^^'-"'^^ F{s, X;9^/'^); 
then G satisfies an equation of the form 
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where the conditions in Setting 2 and the choice of u: are such that A/i, as 
it is seen after straightforward algebra, manifestly preserves analyticity in 
9. Using (IHU, analyticity of t'^^^"^ f{C,,t) in follows provided |C| is large 
enough (depending on T). 

6.5. Proof of Corollary 1441 Substitution gives for fix, t). defined bv l|81|) . 
an equation of the form ©, with m = 1, d = 1. Then in k is scalar. 

The vector q is 3 dimensional, indexed by I = l,j = 1,2,3. The 

nonlinearity is quartic and the equation is linear in the derivatives of /, 
thus the only nonzero values of bq^k are when q is (and k = 1, ...,4) or a 
unit vector ej S (and fc = 0, ...,3). Further, it is found that 

5 2 

Jr. = 1, K = 2, cji = - = /3, 71 = 72 = 1, /3i = 3, /32 = -, n = 3 

and in 1)82(1 we have 
4 

(114) ao,l = aO,2 = -1, "ei,! = 2, aez.l = 1> "eg,! = 

This is sufficient to check that Theorem |31 applies. 

Since large corresponds to |C| = large, and argz G 

(— Ivr, Itt) corresponds to avg( G (— |7r, |7r). Theorem El implies that for 
any (p £ (0' f ) large x £ T>^ and large C = x/t^^^ we have 

\fix,t)\ = O {\xf%N+l)m-^t\x\-'/^) = O (|x|5/3t^+lf5(^+i)(|xr3, |2;|-2/3) 

Changing variables, this implies 

xiz)-'fixiz,t),t)=0(^t''+'\z\-hiN+i)i\z\-'^M~^ 

= o(^t^|z|-H(^)(|2;|~i,|z|-i) 

as needed for asymptoticity. The convergence in the series representation in 
t''/^ follows from Theorem |31 (ii) . It is seen from (|94() that all the exponents 
of t are integer multiples of |. □ 

Note 49. Large includes part of the region where Theorems^ and\^ imply 
Borel summability of the expansion in inverse powers of z. Together, the 
results provide uniform control of the solution. 

7. Appendix 

7.1. Asymptotic behavior: further comments. In the assumptions of 
Theorem 121 by the remark following it, formal series solutions to the initial 
value problem are asymptotic to the actual unique solution. The discussion 
below addresses the issue of deriving this series, or, when less regularity is 
provided and only the first few terms of the expansion exist, how to show 
their asymptoticity. 

Heuristic calculation. Assuming algebraic behavior of f in our assumptions 
on the nonlinearity, it is seen that the most important terms for large x 
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(giving the "dominant balance") are f^, Voi, coming from the constant part 
of V, and r(x, t). This suggests that, to leading order, f(x, t) ~ f/(x) + 
e~^''(*^'^)r(x, r)dr. If we substitute 

(115) f(x,t) = Ai(t)x-"'-i + f 

into 0, f will generally satisfy an equation of the form Q, for an increased 
value of ar] if the process can be iterated, as is the case in the examples in 
|HI, it generates a formal series solution. 

To obtain rigorous estimates, one writes the equation for f defined in 
H115() and applies Theorem^ to show f = o(x~"''-'-). If the coefficients of the 
equation allow it, this procedure can be repeated to obtain more asymptotic 
terms for f. This is the case for instance in the assumptions of Theorem |21 
where a complete series is obtained, which is furthermore Borel summable 
to f. 

The discussion also shows that the assumption Or > 1 can be often be 
circumvented by subtracting the higher powers of x from f . 

7.2. Simple examples of Borel regularization. In this section we dis- 
cuss informally and using rather trivial examples, the regularizing features 
of Borel summation. An excellent account of Ecalle's modern theory of 
generalized summability is found in |13j : see JOl as well. Many interesting 
results, using more classical tools can be found in 'T'. 

Singular perturbations give rise to nonanalytic behavior and divergent 
series. Infinity is an irregular singular point of the ODE f — f = ^/x, and 
the formal power series solution / = X]fet=o(~-'-)^^'-^ '^"^ diverges. In the 
context of PDEs, the solution h of the heat equation ht — hxx = with 
/i(0, x) real-analytic but not entire, has a factorially divergent expansion in 
small t, the recurrence relation for the terms of which is kH^ = H'f^_^. 

The Borel transform of a series, is by definition its term-wise inverse 
Laplace transform, which improves convergence since C~^x~^~^ = p^/kl. If 
the Borel transformed of a series converges to a function which can be con- 
tinued analytically along and is exponentially bounded, then its Laplace 
transform is by definition the Borel sum of the series. Since on a formal level 
Borel summation is CC~^, the identity, it can be shown to be an extended 
isomorphism between series and functions; in particular, the Borel sum of / 
above, C{l+p)~^ is an actual solution of the equation. Another way to view 
this situation is that Borel transform maps singular problems into more reg- 
ular ones. The Borel transform of the ODE discussed is (p + l)£~^/ + l = 0. 
The inverse Laplace transform of ht = h^x in 1/i is hxx — phpp ~ ^hp = 
which becomes regular, Uxx — Uzz = by taking h(p,x) = p~^/^n(2p^/-^, x), 
z = 2pV2. 

It is in its latter role, of a regularizing tool, that we use Borel summation 
in PDEs. 
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7.3. Derivation of equation from (jl]). We define an m-dimensional 
vector f by ordering the set |c^u : < |j| < n|. It is convenient to intro- 
duce g2(x, t, f) so that 

g2,J (x, t, {^U}|j|<„_i) = - ^ g2,i(x, t, {)d:rj 

J|=n i 

So, for showing that @ implies © it is enough to show that for I < n' < n, 
for U'l = n' — 1, 



di' 



gl(x,t,f)+^g2,(x,t,f)a,T 



is of the form on the right hand side of ©• We do so in three steps. 
Lemma 50. Consider for k >1, 



(116) E(x, t) = Y\{^, t, f ) n 

qbO {m;k} 

where {m;k} denotes the set {(/,j) : 1 < / < 1 < |j| < k}, and f means 
summation over q with the restriction 

(117) Yl ^ ^ 

{m;k} 

Then, for i = 1, 2.., d, 9a;.E(x, t) has the same form as Ulb]) with restriction 
jll'^ , provided k is replaced by k + 1. 

Proof. The proof is straightforward, keeping track of the number of deriva- 
tives and the powers involved: note that 



^— bq(x,t,f)a,Jz+a,,bq(x,t,f) n (4/;^'"' 
1=1 / {m;k} 



m k 

+5^bq(x,t,f))^ ^ q,,y [cif.y''' d^Mfi') n i^if^y"' 

qbO l'=l |j'|=l {m;k} 

where indicates that the term I = = j' is missing from the product. 
Manifestly, this is of the form ()116|) with a suitable redefinition of bq and 
with the product of the number of derivatives times the power totaling at 
most 

+ 1 + \j'\{qi',y - 1) + ^IJkij = 1 + Ijkij < ^ + 1 

{m;k} {m;k} 

Hence restriction ()117|) holds, now with k + 1 instead of /c. □ 
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Lemma 51. For any n' > 1, and any J' with |J'| = n' — 1, 
(118) a^'gi(y,t,f(y,t)) = ^^bq(x,t,f) n {^^f\ 

qbO {m;n'-l} 

for some bq, depending on n' , gi, and its first n' — 1 derivatives, and where 
means the sum over q with the further restriction 

YMWi < n' - 1 

{m;n' — 1} 

Proof. The proof is by induction. We have, with obvious notation, 

^a;igl(x,t,f(x,t)) = gi,^^ +gi,f • 

which is of the form ()118|) . Assume H118|) holds for n' = k > 1, i.e. for all 
J' satisfying |J'| = — 1, 

aJ'gi(x,t,f) = ^V(x,t,f) n (^x//)"' 

qhO {m;k-l} 

Taking a xi derivative, and applying LemmaEIl d^Siiv^ *i f ) for | J| = k will 
have the form above, with k — 1 replaced by k and with restriction 

{m;k} 

Thus, H118|l holds for n' = k + 1, with a different b. The induction step is 
proved. □ 

Lemma 52. For n' = 1, 2, n, anc? any J mi/i | J| = n' — 1 we have 



(119) 9^[g2,i'(x,t,f)5x,f] =5]*bq(x,t,f) J] (^dlfi"" 

qhO {m;n'} 

/or some bq, depending on n' , g2 and its first n' — 1 derivatives, where X^q^o 
denotes summation with the restriction 

(120) 

{m;n'} 

Proof. Clearly ()119() with restriction ()12U() holds for n' = 1. Suppose it holds 
for n' = k. Then we note that if |J| = A; + 1, then there exists some index 
\ < i < d and some J', with |J'| = A; so that = dxi[d^ ]; hence applying 
Lemma Em we obtain (fTTU)) and ({1^ for n' ={k + l). □ 

7.4. Some useful inequalities. 

(1) We start with a simple inequality for a > 1 and ^ > 0: 



(121) (1 + /i") /\"-ie-^"ds < 2r(a) 

Jo 

This is clear for /i < 1, while for /x > 1 we write (1 + /x") < 2/i" and 
note that s^-ie-'^^ds = /x-"r(a). 
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(2) For Q > 0, /J, > 0, fj = 0, 1, > 2 and m G N, 

(122) + - + + ^ ^"^ 

where C{m) is independent of a and v. Indeed, the integral is 
bounded by 

\Jo J I /[l + /i2(l-s)2]- - (1+^2/4)- io 



e 



+ max ^— -tt;— / ^ds < ^^'T/, h 



2r(a)(^i/)-° e-/^^/^ 



.e[i/2,i] [l + //2(l-s)2]'T7g - (l + ^2/4)a a(l + /i2/4)'^ 

2r(Q)(^i/)-" 2°+^r(a)(/ii/)-" (^i/)"e-^'^/2 
- (1 + ^2/4)'^ + + ciT+ 2"+iar(a) 

2r(a)(/ii^)^" 2"+ir(a)(^i/)-" 



(3) For n > 1 the function 



is bounded in M"^, as it can be checked applying Watson's lemma for 
large fi and noting its continuity on [0, oo). Thus, for some constant 
C and 1/ > 1 we have 



1^1 ■ - ■ - C\P\ 



(123) / e^Nr+-lf-^l"d, < ^iM^e-br 
^ ' Jo " 1 + IpI" 

(4) We have |p''| < maxj<rf l^l^*^' < X]i<d ^^'^ thus for some con- 

stant C and all j < m we have 

(124) \V,{-p)\<CY,{l + \pr) 

i 

Also, for some C72 > 0, |Pj(-p)| < ^2 + IkI + 1^^!) =: C2(d + (z) 
and thus, for > C2 + 1 we have, for < /' < n, 

(125) IpI'' rel^.(-p)l(*--)e^(-+i)'?dT< |p|''e«'^+^2*'^ Te^^^^^^^^dr 

JO JO 

< 7^1-^7«„^g(t+l)+C2trf IpI' 1 ~ g"'^ ^ ^3(^) ^^q(t+l)+C2fQ! 

[(Z. - C2)q]''/" '^^^ ^l-^'/n - (^ - C2y'/^ 

7.5. Modified estimates for Lemma I38L ^From (|123jl it follows that for 
a constant C independent of $ we have 

(126) I* * *| < Ce'^(*+i)^»(l*''l+l*''l")||*||^„||*||^„ 
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In particular !B(z^, n, S) is a Banach algebra. For the equivalent of Lemmall6| 
we use the following bounds. 

(127) 

1= /''"'5-ie-(*+i)[lPil"-(bil-^)"]e-'^(*+i)^d5 < t"'^ s^~'e-'^'+'>ds 



< 







r(Q)(t + i)« 

and I<\p,r /'s"-^e-^(*+i)IP^I"[i-(i-)"]ci.<C, 2"r(a)|pir 



Ki + l)lRl 

where we used p22|) for o" = 0. From (|127|) it is clear that 

(128) ||H*F,-||,„ < |||H| * < C[Tia)fc''{u{t + l))~''''\\F\\, 

In Lemma 1211 we get instead 



|F| * IGI 



Very similar changes are made in in Lemma OH Corollarv l251 and in Lemma OHl 
where in the proof we use H125|) instead of H45jl . Definition I27L Lemma 1281 
and Definition EHl do not change. Lemma IHHl Lemma IHTl change in the same 
way as above. In Lemma EH we use again (I125p instead of (|45p to make 
corresponding changes. Finally, in Lemma EHl changes to v/A/c. 
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